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1. Introduction 



B meson decays into a pair of light mesons are mediated by rare flavour-changing b — > u 
("tree") or loop-induced b — > d,s ("penguin") quark transitions. The plethora of hadronic two- 
body final states, consisting of e.g. n,p,K,K*,ri or <p mesons, opens a particularly rich laboratory 
for testing the CKM paradigm of flavour mixing and CP violation. 

In order to exploit the wealth of data that has been collected at the B physics experiments, a 
quantitative control of the strong-interaction effects is essential. QCD factorization (QCDF) [JTJ is a 
systematic framework to compute the hadronic matrix elements from first principles. It is based on 
the statement that the matrix elements of the operators in the weak effective Hamiltonian factorize 
in the heavy quark limit mt, Aqcd according to 

{MiM 2 \Qi\B) ~ F BMl (0) f M2 J du T({u) <j) Ml (u) 

+ h Sm x /m 2 J dcodvdu 7/ 7 ((0,v,m) <\> b {(Q) 0m, (v) 0m 2 ("), (1-1) 

where the non-perturbative effects are confined to some process-independent hadronic parameters 
such as decay constants /m, light-cone distribution amplitudes <pM and a heavy-to-light form factor 
F BM (0) at maximum recoil. The short-distance hard-scattering kernels T?' 11 , on the other hand, are 
perturbatively calculable and currently being worked out to next-to-next-to-leading order (NNLO), 
i.e. at ff(ccf) [Q, |L 0, §]. Here we report on the status of the perturbative calculation, which we 
divide into two parts: vertex corrections (7/) and spectator scattering (T? 1 ). 



2. Spectator scattering 

We start with the class of short-distance interactions that involves the spectator quark of the 
decaying B meson. Technically, the calculation of the kernels 7} amounts at the considered order 
q s 2 to a 1-loop calculation with six external legs, cf. the left diagram of Figure [j] The calculation is 
complicated by the fact that the interactions between the soft spectator and the energetic (collinear) 
particles in the final state induce a new (virtual) degree of freedom. These so-called hard-collinear 
or jet modes describe configurations of energetic massless particles with virtualities }ij lc ~ wi^Aqcd, 
which is in between the hard scale }Xh — and the hadronic scale Aq CD . 

The decomposition in ( |L~I| ) relies on a perturbative treatment of these hard-collinear modes 
(it is usually assumed that jj,/ lc ~ 1.5 GeV). As the kernels T/ 1 contain the effects from two short- 
distance modes with different virtualities, the calculation becomes most transparent when it is 
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Figure 1: Two-step matching procedure QCD — > SCETi(hc,c,s) — > SCETn(c,s) of the spectator scattering 
conttibution. S (s) denotes a soft heavy (light) quark, c (c) a collinear quark in the direction of the meson Mi 
(M2) and h (he) refers to virtual hard (hard-collinear) modes. 
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organized as a two-step matching calculation between QCD and soft-collinear effective theory 
(SCET) [g], where the perturbative degrees of freedom are subsequently integrated out. 

The matching calculation is illustrated in Figure fft QCD is first matched at jU ~ jit/, onto an 
effective theory called SCETi(hc,c,s) and then at a lower scale /I ~ /i/ !c . onto SCETn(c,s). The 
respective matching coefficients Hi (J) encode the effects of the hard (hard-collinear) modes, while 
the hadronic matrix element of the remnant non-local operator in SCETn(c,s) yields the factorized 
light-cone distribution amplitudes 0m of the three mesons. The hard-scattering kernels Tf 1 finally 
follow from a convolution of the perturbative coefficient functions 

Ti n {(o,v,u) = J dzJ{(0,v,z)Hi(z,u). (2.1) 

From Figure [j] it is evident that the spectator scattering mechanism requires at least one perturbative 
gluon exchange with J = &{a s ) and Hi = &{\) (at tree level). It has been pointed out in [0] that the 
same jet function / also enters the factorization formula for heavy-to-light form factors [§]. As the 
0(ccj) corrections to / have already been worked out in this context the NNLO calculation of 
the kernels T' 1 reduces to the computation of the &((X S ) terms of the hard coefficient functions 
This programme has recently been completed: the corrections for the topological tree amplitudes 
have been computed in ^ and the ones for the QCD and electroweak penguin amplitudes in |3[]. 

The work in [Q] follows an alternative approach to compute the kernels Tf 1 (for the tree am- 
plitudes). In this work the calculation has been performed in pure QCD by expanding the 1-loop 
diagrams to leading power in 1/m^. The calculation thus yields directly the convolution from 



( |2.1| ) without disentangling hard and hard-collinear effects. The findings of [|4p agree with those 
from [Q, |9|], which demonstrates the equivalence of the diagrammatical approach (QCDF) and the 
effective field theory formulation (SCET) at a non-trivial fixed order in perturbation theory. 



3. Vertex corrections 

The calculation of the kernels T- is conceptually simpler, since the complicated interplay of 
soft and collinear dynamics, which induces the hard-collinear degrees of freedom, is absent in this 
case. The calculation is, however, technically demanding as it amounts to a 2-loop calculation with 
four external legs, cf. Figure ||. 

The 2-loop corrections to the kernels T- for the topological tree amplitudes have been com- 
puted in The work makes use of a couple of advanced techniques, which are widely applied 
in multi-loop calculations. First and foremost the calculation is based on an automatized reduction 
algorithm, which uses integration-by-parts techniques [[K]] to relate the entire 2-loop calculation 
to an irreducible set of scalar Master Integrals (Mis). The actual number of 36 Mis turns out to 
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Figure 2: Matching of the 2-loop vertex corrections (the spectator quark is irrelevant for this contribution 
and not drawn). The different configurations are described in the caption of Figure 1. 
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be large and the presence of soft and collinear infrared divergences in intermediate steps of the 
calculation implies that the Mis are required to up to five coefficients in the £-expansion, which 
makes the calculation somewhat involved (e = (4 — d) /2 in dimensional regularization). 

The whole set of Mis has been computed in Several sophisticated techniques, among 



these the method of differential equations [|1 1|], the formalism of Harmonic Polylogarifhms [ |12| ] and 
Mellin-Barnes techniques [13], have been applied to derive analytical expressions for all Mis. The 
results have further been checked numerically with the method of sector decomposition [14]. About 
half of these Mis have later been confirmed and applied in a calculation for inclusive semileptonic 
B —* X u £v decays by various groups [p3[[. One particular coefficient of one of the most complicated 
Mis has first been given numerically in [|5[ [T3|] and later been rederived in an analytical form [pi]] . 

Apart from the calculation of 2-loop integrals, the NNLO calculation of the kernels T t reveals 
further technical difficulties. The factorization formula ( |1.1| ) implies in particular that the (infrared 
divergent) contributions, which belong to the formal expansion of the non-perturbative objects 
(F BM \/m, <Pm), have to be subtracted from the partonic calculation. In dimensional regularization 
this matching calculation involves non-physical (evanescent) operators. In order to assure that 
these operators disappear from the final factorization formula, their renormalized matrix elements 
have to be made to vanish. The specific pattern of required infrared subtractions is particularly 
complicated in the calculation of the colour-suppressed tree amplitude, which involves a Fierz- 
evanescent operator already at tree level. 

The calculation itself provides a couple of stringent cross-checks. Whereas individual 2-loop 
diagrams contain up to l/e 4 (soft and collinear) infrared divergences, the kernels T* turn out (as 
predicted) to be free of any singularities, which follows after an intricate subtraction procedure 
of ultraviolet and infrared divergences. Moreover, somewhat subtle renormalization scheme de- 
pendent contributions enter at the 2-loop level scale-dependent terms, which can also be checked. 
Another important check of the scheme independence of the NNLO results would be possible given 
a Fierz-symmetric renormalization scheme of evanescent operators. Unfortunately, such a scheme 
has not yet been worked out at NNLO and therefore this final check is currently lacking. 

In view of the complexity of the considered calculation, it is desirable that there are (at least) 
two independent calculations of the kernels T- (as for the spectator scattering kernels T? 1 ). One on- 
going calculation of the tree amplitudes is close to finish [jn|]. The only missing ingredient of the 
NNLO calculation then consists in the calculation of the kernels T/ for the penguin amplitudes [JH 



4. Compilation of NNLO results 
4.1 Conceptual issues 

First of all it is worth noting that factorization has been found in all available NNLO con- 
tributions to work technically, i.e. soft and collinear infrared singularities factorize as predicted at 
NNLO and the resulting convolution integrals are finite (in particular free of endpoint singularities). 

Another important conceptual point deals with the question if the dynamical hard-collinear 
scale H/ 1C ~ y 7 m^AQCD should be treated as a perturbative scale, an issue that has been controversial 
in the literature [J7|, [19[ ]. The question is related to the perturbative expansion of the spectator 
scattering mechanism discussed in Section ||. As this contribution starts at £?(a s ), the NNLO terms 



4 



NNLO corrections to charmless hadronic B decays 



Guido Bell 



constitute the first radiative corrections to this mechanism. The factorization of this contribution 
within SCET, as outlined in Section ||], further allows to systematically disentangle the effects from 
the scales jU/, and H/ 1C and in addition to resum logarithms In \ih /n/ ic via renormalization group 
techniques. The explicit NNLO results [||, |3| £|, ||] do not show any sign of abnormally large 
corrections and the remnant dependence on the scales }Xh and \i] lc turns out to be well-behaved. The 
NNLO calculation thus suggests that the expansion in a s {jXh c ) is under control. 

4.2 Tree amplitudes 

One important subset of hard-scattering kernels, which specify the topological tree amplitudes, 
are by now completely determined to NNLO [^.H, H]. It is interesting to illustrate the structure of 
the perturbative expansion at the amplitude level. Using the input parameters specified in [|20|], the 
colour-allowed (a{) and colour-suppressed ((Xq) tree amplitudes in the B — > Tin channels become 

ai{%%) = 1.008 1 v(0) + [0.025 + 0.010/] y(1) + [0.027 + 0.032/] y(2) 
-0.012| s(1) - [0.021+ 0.015/] s{2) -0.014| p 
= 1-013±SS + (+0.027^1)/, 
a 2 (nn) = 0.223 | y(0) - [0.174 + 0.075/] y(1) - [0.032 + 0.051/] y(2) 
+ 0.090 1 s(1) + [0.034 + 0.025/] s(2) + 0.055 \ p 
= 0.195l»;»* + (-0.101^)/, (4.1) 

where the various terms of the perturbative expansion have been denoted by V^°> corresponding to 
the tree level contribution, (1-loop vertex corrections) and (tree level spectator scattering) 
to the NLO corrections and (2-loop vertex corrections) and (1-loop spectator scattering) 
to the new NNLO terms, whereas the numbers denoted by P give an estimate of power corrections 
~ l/mb that are related to subleading twist wave functions of the pions. 

From < J4- 1| ) it is obvious that the NNLO corrections are particularly important for the imaginary 
parts of the amplitudes and hence for strong phases and direct CP asymmetries, which are first 
generated at 0{oc s ). The new corrections are in some cases found to exceed the NLO terms, which 
can be explained by a numerical enhancement from the Wilson coefficients. In absolute terms, 
however, the new corrections are small and perturbation theory seems to be well-behaved. 

The uncertainties of the NNLO prediction for the colour-allowed tree amplitude (X\ are at a 
satisfactory level of few percent. The colour-suppressed amplitude a 2 , on the other hand, suffers 
from substantial uncertainties, which can be traced back to the strong cancellation between the 
terms denoted by V^°\ and V^. This makes the real part of a 2 particularly sensitive to 
the spectator scattering mechanism, which is proportional to the hadronic ratio fnfB/^Bpf 71 ^)- 
The poor knowledge of the B meson parameter 1/Ag = J °° dco/co #b(<b) in particular makes the 
theoretical prediction of a 2 rather uncertain, which calls for further theoretical progress from non- 
perturbative methods (the numbers are given for Xb = (400 ± 150)MeV). 

4.3 The approximate tree decays B — > K~K°/p~p 

Isospin symmetry implies that the decays B — > n~ n Q / p~ p° are free of QCD penguin con- 
tributions (they depend, however, on small electroweak penguin amplitudes). As they do not re- 
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ceive contributions from weak annihilation either, which constitutes an important class of non- 
factorizable power corrections to the factorization formula ([ET|), their branching ratios are particu- 
larly suited to test the strong interaction dynamics of the tree amplitudes. Treating the electroweak 
penguin amplitudes in the NLO approximation [[]], ^lj ], the QCDF prediction becomes [ ^0| ] 

10 6 Br(B- -> *-*<>) = 6.221};^^ ± o.i6 ± o.43 (5.59j£"), 

10 6 Br(B- -p f p°) = 21.0^ : 9+7.4+o.5+j.5 (22.5±i|), (4.2) 

which is in good agreement with the experimental values given in parentheses. Here L refers to 
the longitudinal polarization and the uncertainties are, in order, due to CKM parameters, hadronic 
parameters, higher order perturbative corrections and power corrections. The theoretical prediction 
depends, however, strongly on the input values for \V u b\ and the form factors Ff n (0) and Aq P (0). 

The theoretical prediction can be made independent of these input parameters by normalizing 
the hadronic decay rates to the differential semileptonic decay rates at maximum recoil. The ratio 

UB ~ K n " } - ~ 3n 2 f 2 \V lld \ 2 \a 1 (nn) + a 2 (nn)\ 2 (4.3) 



dT(8° -^K+i-Vi)/dq 2 



q 2 =0 



provides an exceptionally clean probe of the QCD dynamics of the tree amplitudes [g2J. The 
NNLO prediction for this ratio (0.70+o ^)GeV 2 compares again well to the experimental value 
(O.Sl^Q'J^GeV 2 , which strongly supports the factorization assumption [20]. It would, however, be 
interesting to see if the tendency between the experimental and the theoretical value is reproduced 
in the p -sector. As the semileptonic B — > p£v decay spectrum has not yet been measured precisely, 
one may instead consider a ratio of two hadronic decay rates 

r ( g -> PlPI) „ I«i(PlPl) + «2(PlPl)| 2 

r(SP^p+p£) 2|a 1 (p L p L )P ' K ' } 



which, in contrast to ( |4.3| ), receives corrections from QCD penguin amplitudes and weak anni- 
hilation. The NNLO prediction for this ratio (0.65+q'Jj) is again found to be somewhat smaller 
than the experimental value (0.89^q J^), which may be considered as a hint at somewhat enhanced 
colour-suppressed amplitudes that could be realized in QCDF by a smaller value of the B meson 
parameter Xb — 250MeV 



5. Conclusions 



The NNLO calculation for charmless hadronic B meson decays is particularly important in 
respect of direct CP asymmetries that are first generated at ff{a s ). It opens in particular a new 
mechanism from spectator scattering that induces strong phases and settles some conceptual aspects 
that bring the factorization framework onto a more rigorous footing. 

Whereas the topological tree amplitudes are by now completely determined to NNLO, the 
computation of the penguin amplitudes is to date still incomplete. Further improvements on the 
calculation require a better understanding of power corrections, in particular on the role of the phe- 
nomenologically important scalar penguin amplitude, and more precise determinations of hadronic 
input parameters. 
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